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Abstract 

We introduce a new class of self-similar Gaussian stochastic processes, where the covariance is 
defined in terms of a fractional Brownian motion and another Gaussian process. A special case is 
the solution in time to the fractional-colored stochastic heat equation described in Tudor (2013). 
We prove that the process can be decomposed into a fractional Brownian motion (with a different 
parameter than the one that defines the covariance), and a Gaussian process first described in Lei 
and Nualart (2008). The component processes can be expressed as stochastic integrals with respect 
to the Brownian sheet. We then prove a central limit theorem about the Hermite variations of the 
process. 


1 Introduction 


The purpose of this paper is to introduce a new class of Gaussian self-similar stochastic 
processes related to stochastic partial differential equations, and to establish a decom¬ 
position in law and a central limit theorem for the Hermite variations of the increments 
of such processes. 

Consider the d-dimensional stochastic heat equation 


du 1 

-— = -Am + W^ t > 0, X G M'^, (1.1) 

at 2 

with zero initial condition, where ID is a zero mean Gaussian held with a covariance of 
the form 


E 


W^{t,x)W^{s,y) 


7o(t - s)A(x - 1/), s, f > 0, X, 1/e 


We are interested in the the process U = {Ut, t > 0}, where Ut = u{t, 0). 

Suppose that W is white in time, that is, 70 = Sq and the spatial covariance is the 
Riesz kernel, that is, A(x) = Cd,/3|x|“^, with (3 < mm{d, 2) and Cd,^ = /2) /T{{d— 
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(3)/2). Then U has the covariance (see [T5] i 


E[UtUs] = D(^{t + sY-^-\t-s\^-^y s,t>0, (1.2) 

for some constant 

(1.3) 

Up to a constant, the covariance fll.2p is the covariance of the bifractional Brownian 
motion with parameters H = ^ and K = 1 — We recall that, given constants 
H E (0,1) and K E (0,1), the bifractional Brownian motion = {Bl^’^,t > 0}, 

introdnced in [3], is a centered Gaussian process with covariance 

Rh,k{s, t) = ^ -\t- , s, t > 0. 

When it' = 1, the process B^ = B^'^ is simply the fractional Brownian motion (fBm) 
with Hurst parameter H E (0,1), with covariance RH{s,t) = RH,i{s,t). In [1], Lei 
and Nualart obtained the following decomposition in law for the bifractional Brownian 
motion 

= C^B^^ + 

where B^^ is a fBm with Hurst parameter HK, the process is given by 




- e-y^)dWy, 


(1.4) 


with W = {Wy,y > 0} a standard Brownian motion independent of B^’^, and Ci,C 2 

are constants given by Ci = 2^~^ and C 2 = sj y{i-k) • process Y^ has trajectories 

which are inhnitely differentiable on (0, cxd) and Holder continuous of order HK — e in 
any interval [0,T] for any e > 0. In particular, this leads to a decomposition in law of 
the process U with covariance fll.2p as the sum of a fractional Brownian motion with 
Hurst parameter | — f plus a regular process. 

The classical one-dimensional space-time white noise can also be considered as an 
extension of the covariance (II.2p if we take (3 = 1. In this case the covariance corresponds, 
up to a constant, to that of a bifractional Brownian motion with parameters H = K = ^. 

The case where the noise term hU is a fractional Brownian motion with Hurst pa¬ 
rameter H E (|, 1) in time and a spatial covariance given by the Riesz kernel, that 


is. 


E 


W^{t,x)W^{s,y) 




where 0 < /3 < mm{d, 2) and an = H{2H — 1), has been considered by Tudor and Xiao 
in [13]. In this case the corresponding process U has the covariance 




u — V 


‘^{t + s — u — v) '^dudv. 


E[UtUs] = Dan 


(1.5) 
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where D is given in fll.Sp and 7 = This process is self-similar with parameter H — ^ 
and it has been stndied in a series of papers [D m [m Ha E3]. In particnlar, in |T3] 
it is proved that the process U can be decomposed into the snm of a scaled fBm with 
parameter H — and a Gaussian process V with continuously differentiable trajectories. 
This decomposition is based on the stochastic heat equation. As a consequence, one can 
derive the exact uniform and local moduli of continuity and Chung-type laws of the 
iterated logarithm for this process. In [n], assuming that d = 1, 2 or 3, a central limit 
theorem is obtained for the renormalized quadratic variation 

n—1 

Vn = {(f/(,+i)T/n - U^T/n? " E " U.T/nf] } , 

j=0 


assuming | < TT < |, extending well-known results for fBm (see for example [5l Theorem 
7.4.1]). 

The purpose of this paper is to establish a decomposition in law, similar to that ob¬ 
tained by Lei and Nualart in |1] for the bifractional Brownian motion, and a central limit 
theorem for the Hermite variations of the increments, for a class of self-similar processes 
that includes the covariance dESp. Consider a centered Gaussian process {Xt,t > 0} 
with covariance 


R{s,t)=E[XsXt\ 


E 


Zt-rdBl^ 


Zg-rdB, 


H 


( 1 , 6 ) 


where 


(i) B^ = {B^,t > 0} is a fBm with Hurst parameter H e (0,1). 

(ii) Z = {ZtR > 0} is a zero-mean Gaussian process, independent of B^, with covari¬ 
ance 

E[Z,Zt] = {s + t)-\ (1.7) 

where 0 < 7 < 2H. 

In other words, X is a Gaussian process with the same covariance as the process 
{/q Zt-rdB^R > 0}, which is not Gaussian. 

When H G (|,1), the covariance (11.61) coincides with (II.5p with D = 1. However, 
we allow the range of parameters 0 < id < 1 and 0 < 7 < 2H. In other words, up 
to a constant, X has the law of the solution in time of the stochastic heat equation 
(EHi. when H G (0,1) and d > 1 and (3 = d — 27. Also of interest is that X can be 
constructed as a sum of stochastic integrals with respect to the Brownian sheet (see the 
proof of Theorem 1). 


1.1 Decomposition of the process X 

Our hrst result is the following decomposition in law of the process X as the sum of a 
fractional Brownian motion with Hurst parameter | = id — ^ plus a process with regular 
trajectories. 
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Theorem 1. The process X has the same law as +Yt^t> 0}, where 

dz, (1.8) 


1 r 


K = 


r(7)P 1 + ^2 


is a fBm with Hurst parameter Oij2, and Y (up to a constant) has the same law as 
the process Y^ defined in (El, with K = 2a + 1, that is, Y is a centered Gaussian 
process with covariance given by 


E [YtY,] = Ai / dy, 

Jo 


where 


Ai = 


An 


fOO ^l-2H 


T{-f)T{2H + l)sm{nH) Jo 1 + 
The proof of this theorem is given in Section 3. 


dp. 


1.2 Hermite variations of the process 

For each integer q > 0, the gth Hermite polynomial is given by 

^2 d'^ 

HJx) = 

^ ^ ^ dxi 

See [HI Section 1.4] for a discussion of properties of these polynomials. In particular, it is 
well known that the family {^Hg, g > 0 } constitutes an orthonormal basis of the space 
L^(M, 7 ), where 7 is the iV(0,1) measure. 

Suppose {Zn,n > 1} is a stationary, Gaussian sequence, where each follows the 
A^(0,1) distribution with covariance function p{k) = E[Z„Z„+fc]. If 
it is well known that as n tends to inhnity, the Hermite variation 


K, = 


E 


(1,9) 


1=1 


converges in distribution to a Gaussian random variable with mean zero and variance 
given by This result was proved by Brener and Major in [2]. In 

particular, if B^ is a fBm, then the sequence {Zj^ni 0 < j < n — 1} defined by 




B 


H 


is a stationary sequence with unit variance. As a consequence, i! < 1 — we have that 


n—1 


= Hf+i -Bf 

n V V ~ n 

j=0 


1 
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converges to a normal law with variance given by 

aj = 1^ ^ (|m + - 2|mp^ + |m - . (1.10) 

m^’L 


See [2] and Theorem 7.4.1 of [5]. 

The above Brener-Major theorem can not be applied to our process because X is not 
necessarily stationary. However, we have a comparable result. 


Theorem 2. Let q > 2 be an integer and fix a real T > 0. Suppose that a < 2 — - . For 


t G [0, T], define, 


[nt\ — l 

Fn{t) = n~^ ^ Hq 

j=0 


( 


AX, 




AX, 


L2(0)/ 


where Hq{x) denotes the qth Hermite polynomial. Then as n ^ oo, the stochastic 
process {Fn{t),t G [0,T]} converges in law in the Skorohod space Z1([0,T]), to a scaled 
Brownian motion {aBt,t G [0,T]}, where {Bt,t G [0,T]} is a standard Brownian motion 
and a = y/o^ is given by 


a 


2 




— 2|m|“ + \m — 1|")^ . 


( 1 . 11 ) 


The proof of this theorem is given in Section 4. 


2 Preliminaries 

2.1 Analysis on the Wiener space 

The reader may refer to [5l[6] for a detailed coverage of this topic. Let Z = {Z(h), h E TL] 
be an isonormal Gaussian process on a probability space indexed by a real 

separable Hilbert space TL. This means that Z is a family of Gaussian random variables 
such that E[Z(/i)] = 0 and K[Z{h)Z{g)] = {h,g)q^ for all h,g eFL. 

For integers g > 1, let denote the gth tensor product of Ti, and denote the 
subspace of symmetric elements of 

Let {e„,n > 1} be a complete orthormal system in Ft. For elements f,gE and 
p G {0,..., q], we dehne the pth-order contraction of / and g as that element of 
given by 

OO 

f®pg= ^ ( 2 . 1 ) 

i\ ,...,2p = l 

where f ®q g = f Zt g. Note that, if /, g G Fi®^, then f ®qg = (/, g)-^Qg- In particular, if 
f,g are real-valued functions in = L^(M^,H^,/i) for a non-atomic measure p, then 
we have 

f®ig= [ f{s,ti)g{s,t2) ^iids). (2.2) 

Jr 
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Let 1-iq be the gth Wiener chaos of Z, that is, the closed linear subspace of 
generated by the random variables {Hq{Z{h)), h ||h||^ = 1}, where Hq{x) is the gth 
Hermite polynomial. It can be shown (see [H Proposition 2.2.1]) that if Z,Y iV(0,l) 
are jointly Gaussian, then 




p! (E [ZY]Y if p^q 
0 otherwise 


For g > 1, it is known that the map 


I,{hn = H,(Z(h)) 


(2,3) 


(2.4) 


provides a linear isometry between (equipped with the modihed norm ^/^.\\ ■ 
and Hq, where Iq{-) is the generalized Wiener-Ito stochastic integral (see [HI Theorem 
2.7.7]). By convention, 770 = 1^ and Io{x) = x. 

We use the following integral multiplication theorem from [6l Proposition 1.1.3]. Sup¬ 
pose / G 77®^ and g G 77®*^. Then 


p/\q 


r’=0 


fp+q—2r (/ ®rg^i 


(2.5) 


where f®rg denotes the symmetrization of / Gr g. For a product of more than two 
integrals, see Peccati and Taqqu [8]. 


2.2 Stochastic integration and fBm 

We refer to the ‘time domain’ and ‘spectral domain’ representations of fBm. The reader 
may refer to PCD] for details. Let £ denote the set of real-valued step functions on M. 
Let denote fBm with Hurst parameter H. For this case, we view as an isonormal 
Gaussian process on the Hilbert space ij, which is the closure of £ with respect to the 
inner product {f,g)s^ = E [/(/)/(g)]. Gonsider also the inner product space 

Ah = |/ : / e n(R), ^ |.F/({)n{|‘-“d^ < ooj , 

where Bf = f{x)e^^^dx is the Fourier transform, and the inner product of Kh is given 
by 

(/. 9 ) a , = -Y [ ( 2 . 6 ) 

Jr 

1 

where Ch = ( r( 2 B+i)lin(qrB) ) ' ■ is known (see P Theorem 3.1]) that the space Ab is 

isometric to a subspace of Sj, and Ab contains as a dense subset. This inner product 
fl2.6p is known as the ‘spectral measure’ of fBm. In the case H G (|, 1), there is another 
isometry from the space 

r poo poo A 

|A//| = |/; j J \f{u)\\f{v)\\u-vf^~‘^dudv<oo> 
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to a subspace of where the inner product is defined as 


{/. 9 )|A,l=-ff( 2 i/-l) 


f{u)g{v)\u — vf^ ^du dv, 


'0 JO 


see [9] or [HI Section 5.1]. 

3 Proof of Theorem [1] 

For any 7 > 0 and A > 0, we can write 


1 r°° 




where F is the Gamma function defined by F( 7 ) = y'^ ^dy. As a consequence, 

the covariance fll.7l) can be written as 


^[ZgZt] — 


r(7) 


yy-l^-(t+s)y^y_ 


(3.1) 


Notice that this representation implies the covariance fll.7p is positive dehnite. Taking 
first the expectation with respect to the process Z, and using formula fl3.ip . we obtain 




I ey^dB^j i^j e^^dB^ 


yj-l^-(t+s)y^y 


1 

mJo 


(e^n[0,i](n), eyn[0,s]{v)),^ 

Using the isometry between Ah and a subspace of (see section 2.2), we can write 


(e^n[o,](n),e^’^l[o,.](n)>. = / |er-^^(-Fl[o,]e^ )(-Fl[o,.]eJ^-) 




^ ^-2 f 


where (J^lp^tje"^') denotes the Fourier transform and Ch = y r( 2 H+p 8111 ( 71 ^) ^" • This 
allows us to write, making the change of variable ^ = yy, 


R{s,t) = 


mcj, 


y 


H JO • 

1 


mcl 


y 


H JO 


y^ + e 

1 + 7^ 


- e-y^) d^ dy 

y^jvvt _ ^-yt^ _ ^-ys^ (3 2 ) 


where a = 2H — 7 . By Euler’s identity, adding and subtracting 1 to compensate the 
singularity of y~°‘~^ at the origin, we can write 


Jyyt _ e-y^ — 


e = {cos{riyt) — 1 + i sm{T]yt)) + (1 — e ^*). 


(3,3) 





















Substituting fl3.3p into fl3.2p and taking into account that the integral of the imaginary 
part vanishes because it is an odd function, we obtain 

2 roo f-oo „l-2ir , 

+ sm{riyt) sm^rjys) + {cos{riys) — 1)(1 — e“^*) + {cos^rjyt) — 1)(1 — e~^^) 

+ (l-e-^*)(l-e-*'")) dydy. 

Let = {B^^\r],t),r] > 0,t > 0}, j = 1,2 denote two independent Brownian 
sheets. That is, for j = 1,2, B^^"> is a continuous Gaussian held with mean zero and 
covariance given by 

= min(r7, x min(t, s). 

We dehne the the following stochastic processes: 


Ut = 


14 = 


Y,= 


^/2 


(*00 poo 


y/T{'y)CH Jo Jo 
V2 


_^_1^ 

y 2 2, 


oo poo 


^/T^ChJo Jo 
V2 

Vmc 


_^_1^ 

y 2 2, 


^l-2H 


1 + 1]“^ 


^1-2H 


1 + 1]“^ 


{cos{yyt) - 1) B^^\dy,dy), 


{sm{yyt))B^^\dy,dy), 


OO poo 


_^_1^ 

y 2 2, 


.1-2H 


H JO JO 


l + y" 




(3.4) 

(3.5) 

(3.6) 


where the integrals are Wiener-Ito integrals with respect to the Brownian sheet. We 
then dehne the stochastic process X = {WG > 0} by W = Gt + 14 + 14, and we have 
E [XsW] = R{s,t) as given in fl3.2p . These processes have the following properties: 

(I) The process 144 = f/j + 14 is a fractional Brownian motion with Hurst parameter | 
scaled with the constant In fact, the covariance of this process is 


E[mW,] = 


oo /•oo „1—2H 

/ 9-“-'' 


Jo Jo I + rj 

+ sm{r]yt) sm{r]ys)^dr]dy 

1 r r 


^(cos(?7|/t) — l){cos{riys) — 1) 


r(7)4^^ 


2 /' 


Mi _ 


H JO 


y^ + e 


l)(e 


—i^s _ 


l)d^dy. 


Integrating in the variable y we hnally obtain 


E[W/W,] 



- l)(e-*^" 


1 ) 




where Ci = j^dz = Taking into account the Fourier transform represen¬ 

tation of fBm (see m page 328]), this implies k, zfF is a fractional Brownian motion 
with Hurst parameter 
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(II) The process V coincides, up to a constant, with the process introduced in fll.41) 
with K = 2a -h 1. In fact, the covariance of this process is given by 


E[y)n] = 


2C2 


r(7)ci 




(3,7) 


H JO 


where 


poo 

C2= 7 -— i;d7]. 

Jo 1 + 


Notice that the process X is self-similar with exponent This concludes the proof 
of Theorem [TJ 

4 Proof of Theorem [2] 

Along the proof, the symbol C denotes a generic, positive constant, which may change 
from line to line. The value of C will depend on parameters of the process and on T, 
but not on the increment width n~^. 

For integers n > 1, dehne a partition of [0, cxd) composed of the intervals {[^, > 

0}. For the process X and related processes U, V, W, Y dehned in Section 3, we introduce 
the notation 

AXi = Xi+i - Xi and AXq = Xi, 

n n n n 

with corresponding notation for U, V, W, Y. We start the proof of Theorem [2] with two 
technical results about the components of the increments. 


4.1 Preliminary Lemmas 

Lemma 3. Using above notation with integers n>2 and j, k > 0, we have 


AWj AhFfe 


n n 


(a) E 
in / li.gl) . 

(b) For j + k>l, 


= In " (|j — fc — 1|" — 2\i — + \j — k — 1|"), where n is defined 


E 


AYi AYk 


n n 


< Cn-^^ij + k) 


a-2 


for a constant C > 0 that is independent of j, k and n. 

Proof. Property (a) is well-known for fractional Brownian motion. For (b), we have from 

dXTI): 


E 


AYj AYk 


j_LAl k 

n n 


2 Co 


r(7)C|,n“ 

2c: 

fioFi 


1 rl 


Jo 


y 


—a+l 



^-y^j+k+u+v)^^ dn dy. 


0 Jo 
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Note that the above integral is nonnegative, and we can bonnd this with 


E 


AYi AYk 


n n 


< Cn 


y-a+l^-y{j+k) 


= Cn-^{j + 

Jo 

< Cn-^{j + ky-^. 


□ 


Lemma 4. For n> 2 fixed and integers j,k > 1, 


E 


AWj AYk 


n n 


< Cn-^f^-^k-'^ 


for a constant C > 0 that is independent of j, k and n. 

Proof. From fl3.4p - (3.6) in the proof of Theorem (H observe that 


E 


AhF, AYk 


n n 


= E 


{AU^+AVyAYk 


n n 


= E 


AUj AYk 


_ 

n n 


Assume s, f > 0. By self-similarity we can dehne the covariance function -0 by E = 

s"E [Lf/sTi] = s“0(f/s), where, using the change-of-variable 9 = px, 


^p{x) = 


r-oo koo „1—2ir 

:y {cos{yr]x) - 1) (1 - e“*^) dp dy 


>0 JO 
poo 

—a—1 


l + T 


y-'^-yi-e-y 


kOO 0l—2H^2H 


'0 


X‘ 


+ 02 


(cos(|/0) — 1) dO dy. 


Then using the fact that 


0l-2Hx2H 


X‘ 


+ 02 


<1^ 


-2H| U|211-l 


we see that |'0(a:)| < Cx'^^ and 

poo 

fj'{x) = 2H y-'^-yi-e-y 


- 2 / y-'-yi-e-y 


fOO q1—2H^2H—1 

) a;2 + 02 

poo Ql—2H^2H+l 


/o (x2 + 02)' 


cos{y9) — 1) d9 dy 
cos{y9) — 1) dO dy. 


Using fld.ip and similarly 


q 1-2H ^2H+l 


(X2 + 02)2 


<1^ 


-2H\ I |2H-2 

\x\ 


we can write 


(4.1) 


(4.2) 


\'f'{x)\<x^^ ^\2H — 2\ / 1/ " ^(1 —e *') / 0 ^ (cos(|/0) — 1) d9 dy < Cx 


^2H-2 
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By continuing the computation, we can find that \'ip''{x) \ < We have for j, /c > 1, 


E 


/\Uj AW 


k 

n n 


= n-'^{k + 1 )“ 

— ('ll! 


J + 1 
k + 1 
J + 1 




T 


k + 1 


= n-“((fc + l)"-r) ( 'ijj 
' 3 + 1 


3 + 1 
k + 1 


-■0 


3 




k + 1 




k+1 


-■0 


k + 1 
J + 1 


+ iIj 


With the above bounds on -0 and its derivatives, the first term is bounded by 




0 


3 + 1 
k+1 

rl 


'ijj 


3 


k + 1 


< an, ^ {k + m)“ ^du 


1 

k+1 


'0 


0 ' 


< Cn-^k' 


k+1 
i-2 f 3 


+ V 


dv 


2H-2 




and 

n-^k^ 




J +1 

k+1 

= n-^k' 

< n-^^k^ 




k + 1 
1 

* k+1 




J + 1 




V-' 


k+1 


+ M I du 


0 ' ( ^ + M ) du 


' 1 
k + 1 


0 ' 


+ u 


du + 


- 1 „i 
k+1 I k 


-7 

fc +1 


< Cn-’^k^-^ 

This concludes the proof of the lemma. 


2H-2 

I ) + Cn-^k^-^3 


| 0 "(m + n)| dv du 

2H-3 


< Cn-^k-^f^-^. 


□ 


4.2 Proof of Theorem [2] 


We will make use of the notation 


have 



AAi 

n 


. From Lemma [3] and Lemma 0] we 

L^{n) 


K,n “(l + 0 j-„). 


where \9j^n\ A if j > 1. Notice that, in the definition of Fn(t), it suffices to 

consider the sum for j > uq for a fixed hq. Then, we can choose Hq in such a way that 
CnQ~^ < i, which implies 

/3l^ > - Cr-^) (4.3) 
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for any j > tiq. 

By (I23D, 







\i. i+i') 

I- n ’ n ' 



5 


where denotes the multiple stochastic integral of order q with respect to the process 
X. Thus, we can write 


lnt\ — 1 

F„{t) = n-i ( Ifl, 

j=no 


L .7 + 1 ' 

n ’ n ^ 


The decomposition X = W + V leads to 


rX I 1 

-'g 1 i±l' 


E 

r=0 


^ I ( -\®r \ I 1 

r/ I '*'[ 1,1+1^ 


We are going to show that the terms with r = 0,...,g — Ido not contribute to the limit. 
Dehne 

GnW=n-= E 

j=no \ n n 


and 

[nt\ — l 

Gr,{t)=n--2 11 ^^^ 
j=no 

Consider the decomposition 


“9 tW I 1 S9C 


Fn{i) — {Fn{t) — Gn{t)) + {Gn{t) — Gn{t)) + Gn{t). 


Notice that all these processes vanish at t = 0. We claim that for any 0 < s < t < T, 
we have 




E[\F^{t) - Gr,{t) 

- (.F„(s) 

-Gn 

(5))P] 


n 

(4.4) 

and 











E[\Gn{t) - Gn{t) 

- (Gn(s) 

-Gn 

(^))r] 

< (M 

— [nsj)*^ 

? 

n 

(4.5) 

where 0 < 

(5 < 

1. By Lemma [3l 

1 |AW,/„|| 

2 

L2(0) 

= nn 

" for every j. As 

a consequence. 

using fj2.4p 

we 

can also write 









[nt\- 

L 







Gn{t) = 

n--2 ^ 



"^awA . 






j=no 







Since k 2 W is a fractional Brownian motion, the Brener-Major theorem implies that the 
process G converges in Zi)([0,T]) to a scaled Brownian motion G [0,T]}, where 
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cr^ is given in fll.lip . By the fact that all the p-norms are equivalent on a hxed Wiener 
chaos, the estimates fl4.4|] and fl4.5p lead to 

E[\F^{t) - G^it) - (Fnis) - Gn(s))p?’] < (NJ (4.6) 

and 

E||G„(i) - g„(() - (G„(i>) - G„(i>))p>’] < (4.7) 

for all p > 1. Letting n tend to inhnity, we deduce from fl4.6p and 04.71) that the sequences 
Fn — Gn and — Gn converge to zero in the topology of Zi)([0, T]), as n tends to inhnity. 

Proof of (fg .- We can write 


5-1 


E [|F„(i) - G„(i) - (F„(s) - G„(s))p] < C^E[«>?,„], 


r=0 


where 


[ntj — 1 




j=[nsJVno 


I (S)r 


tY 


j_ 2+11 j q—r \ [i j + 1 


We have, using 04.3p . 

E|-S>" ] < 


[nt\ — 1 

X E 

j,fc=[nsJVno 


E 


W I 1 0r 


tY 


1 ear I 

■‘■''i .7+1 1 j ^q—r \ '^\i. .7+1 ' 
n’n// \ In’n 


W { -t 


fc+1 \ I -^q—r 

In’ n 


rV 


^q—r 
r_fc k±i' 

L Ti ’ n ^ 


Using a diagram method for the expectation of four stochastic integrals (see m), we hnd 
that, for any j, k, the above expectation consists of a sum of terms of the form 



AW^AWk 

n n 


ai 



AW AW 

n n 



AW^AYk 

n n 


aa 



AW AWfe 

n n 


(I4 




where the a* are nonnegative integers such that ai + 02 + 03 + 04 = g, oi < r < g — 1, 
and 02 < g — r. First, consider the case with 03 = 04 = 0, so that we have the sum 


n 


-l-\-qa 


\nt\ — 1 

E 


E 


AW7 AWfe 


E 


AW AW 


q-ai 


j,fc=[nsJVno 

where 0 < Oi < g — 1. Applying Lemma [3], we can control each of the terms in the above 
sum by 

- A: + 1|“ - 2|j - A;|" + |j - A: - 
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which gives 


n 


-1+qa 


j,k= [nsj Vno 


E 




n n 


ai 


E 


AY, AYk 


j_L^l k 

n n 


q-ai 


< Cn 


-1 


lnt\ — 1 


[nt\ —1 

j,k=lnsiVno,jjtk 


\j= Vn-o 
[nt J — 1 

< Cn-^ ^ (j“-2 + /(“-2)+i) 

j = [nsJ Vno 

< Cn-^ {[nt\ - + {[nt\ - [nsJ)[''(“-2)+2lvo) . 


(4.8) 


Next, we consider the case where 03 + 04 > 1. By Lemma [3], we have that, up to a 
constant C, 


E 


AY, AYk 


3_i-Al _fc 

n n 


< C 


E 


Ahh, Ahh. 


j_L^yv 

n n 


so we may assume 02 = 0, and have to handle the term 


n 


-l+qa 


[nt J — 1 

E 


E 


AW, Ahh. 


j_L^vy 

n n 


q—a3—a4 


E 


Ahh, AY, 


3_i-Al k 
n n 


as 


E 


AF, AWk 


^i-Ayy 
n n 


a4 


j,k= [nsj Vno 

for all allowable values of 03, 04 with 03 + 04 > 1. Consider the decomposition 

[nt J — 1 


n 


— l+qa 


E 


j,k= [nsJ Vno 


= 


E 


AIF, AWk 


n n 


q—a^—a^ 


E 


AIF, AF 


k 

n n 


as 


E 


AF, AWk 


j_/-\yy 

n n 


<24 


[nt J — 1 

E 


E 


AIF? 


1= [nsj Vno 
[ntJ-1 j-l 

E E 

jf=[n5jVno A:=[nsJVno 
[ntJ-1 k-1 


q—a^—a4 


E 


AlFi AF 


J_L^1 J_ 

n n 


+ n' 


qa—l 


E E 

k=[ns\Vno j' = [nsJVno 


E 


E 


AfFi AWk 


^L^yy _fc 
n n 


AW I AWk 

n n 


q—a^—a4 


q—as—a4 


E 


E 


a3+a4 


AW^ AYk 

n n 


AW^AYk 

n n 


as 


as 


E 


E 


(4.9) 


AF,' AWk 


j_L^yy k 

n n 


AY, AWk 


j_L^yy k 

n n 
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We have,by Lemma [3] and Lemma IU 


\nt\ — 1 


n 


-l+qa 


E 


q—az—a^ 


AWiAW 

n n 


as 


E 


AW AWfc 


j_L^yv _fc 

n n 


E AWj AWfc E 

j,k= [nsj Vno 

[nt\ — l 

< Cn-^ E ^■(a3+a4)(a-2) 

j = [nsJVno 

\nt\-l j-1 

a3(2_H'-2)-a47 \ ^ L-a37+“4(2-ff-2) _ i,\{q-a3-a4){a-2) 


\j-k\^ 

_ j|('J-“3-a4)(o-2) 


+ Cn-^ 

ji=[nsJVno /c=[nsJVno 

k-\ 

_|_ Cn~^ ^-a37+a4(2H-2) ja3{2H-2)-a4'y 

k=[ns\\/no j' = [Ti5jVno 

< Cn-i (([ntj - + {[nt\ - [nsJ)[''(“- 2 )+ 2 lvo 

+ ([ntj - L 7 isJ)L 3 ( 211 - 2 )-a 47 +l]V 0 ^ ^[a 4 ( 2 H- 2 )-a 37 +l]V 0 ^ _ 

Then fl4.8p and fId.lOp imply fl4.4p because a < 2 — 


Proof of We have 


[nt\ — l 


L/fcDj—i , 

Gn{t) - Grrit) = n--^ hf 

j=no ^ 


-q 

n 


AW. 


-q 


L^(n) 


jW I -t^q 
^q V ■*■[2- i+i) 


and we can write, using fl4.3p for any j > rio. 


d”" - 

AW. 

-q 

= (k ^n ")2 

(1 + ^j,n) 2—1 

1 JyTL 

n 

L 2 ( 0 ) 






This leads to the estimate 


E 


Gn{t)-Gn{t)-{Gn{s)-Gn{s)) 

\nt\ — 1 


< Gn 


-1 


X 


E r-^+ 


\nt\ — 1 

E 1^ - 

il=[nsJVno j,fc=[nsJVno 

< Gn-^ {[nt\ - + ([ntj - [nsJ)['?(“- 2 )+ 2 lvo) , 


which implies fl4.5p . 

This concludes the proof of Theorem |2l 
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